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We show that the coupled two-fluid gravitating system (e.g. stiff matter and 'vacuum 
energy') could trap nonlinear gravitational waves (e.g. Einstein-Rosen waves). The grav- 
itational wave amplitude varies harmonically in time transferring the energy coherently 
to the stiff matter wave, and then the process goes to the backward direction. This pro- 
cess mimics the behaviour of trapped electromagnetic waves in two-level media. We have 
defined the limits for the frequency of this energy transfer oscillations. 

Keywords: Cylindrical Gravitational Waves; Stiff Matter 

04.30.-w, 47.35.Bb, 98.80.Cq 

Investigation of the nonlinearly coupled state of electromagnetic waves with 
two o r thre e level medium has become one of the hottest research spots this last 
The most exciting feature of this nonlinear effect is that the electro- 
magnetic wave, which has a fixed propagation velocity in vacuum, could form with 
medium some breather-like structures localized in sp ace. S uch systems allows not 
only to slow down, and eventually storj^ a light pulse 5 6 7 , but also there exists a 
possibility to release it, thus creating a gap soliton memory^. In the s imple st case 
of two-level medium the governing equations are Maxwell-Bloch syste m^™ and, as 
the very recent theoretical investigations showf^, one can control the mobility of 
the light pulse by the stationary radiation with a different polarization. 

We want to find out if similar effect exists for nonlinear gravitational waves 
inside of gravitating multi-component fluid systems, which are analogous to multi- 
level mediums in nonlinear optics. 

Existence of gravitational waves is generally accepted, at least since their in- 
direct observation via energy and angular momentum losses in binary pulsars (see 
e gfEl an( j references therein). The study of gravitational waves involves mainly 
the weak field approximation of Einstein's equations (in close analogy with what 
is normally done in electromagnetism), which is obviously the relevant regime for 
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existing gravitational wave detectors. On the other hand, the non linearity of the 
gravitational field is one of its most characteristic properties. It is likely that some 
of the crucial properties of the gravitational field show themselves only through the 
nonlinear terms . Thus , a number of exact nonlinear gravitational wave solutions 
have been found™. 

it was shown that standing gravitational waves do not exist for the main 
classes of exact, wave-like solutions to the Einstein equations. The only exception 
is the Einstein-Rosen solution, which will be considered in this paper. The exact 
nonlinear plane-symmetri c st anding gravitational wave solution, but bounded with 
domain wall, was found irPSl. 

Einstein-Rosen cylindrical waves are amon g the simplest non-stationary solu- 
tions to the vacuum Einstein equation d^-^ i Despite the fact that cylindrically 
symmetric waves cannot describe exactly the radiation fro m bounded sources, they 
have played an important role in many physical models^! Also we note that cylin- 
drical gravitational waves are related to a special class of spherical and toroidal 
wave d 1 ? 1 18 1 . 

One of the important features of the Einstein-Rosen wave is that it is special 
case of the gravitational Bessel beairP22. Bessel beams are cylindrically symmetric 
solutions to the wave equations, which were exte nsively studied and applied in 
both optics and acoustics (see the reviewd^- ^ 1 ^^ . The peculiarity of this type of 
wave is that it is well localized in space. This is the property we are looking for in 
interactions of gravitational waves with media. 

There have been numerous investigations o n the int eraction between multi- 
component plasma and linear gravity waves (see^IMUS] an( j re f erences therein). 
In this paper we consider the nonlinear case. For the simplicity we assume coex- 
istence of two perfect fluids with zero relative velocity that have their own energy 
momentum tensors of the form: 

= (e* + V l Wu v - g^p* ,i = l,2 (1) 

where e l are the energy densities, p l are the pressures and u v is the normalized 
4- velocity vector, 

u"« P = 1 • (2) 

As one of the ingredients of our two-fluid system we take the negative vacuum 
energy, 



e 



P=-P, (3) 



with the energy-momentum tensor 

T% = -W . ( 4 ) 

which can be associated to one or more scalar fieldJ^. Having in mind the elec- 
tromagnetic analogy we assume that under the influence of nonlinear gravitational 
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wave the vacuum energy p can oscillate around it's average value (but stay always 
positive) . 

For the second fluid we denote e 2 = e, p 2 = p and consider the equation of state 

e = p + 2p, (5) 

which for the case of p — > reduces to the equation of state of the stiff matter, 

e = p. (6) 

The study of stiff fluid in gravitational waves contexts is interes ting in the sense that 
in such matter the speed of sound is equal to the speed of lig The equation 

of state ([5]) corresponds to almost maximum stiffness compatible with causality and 
perhaps is good approximation for supernuclear densities of matter existing at early 
stages of the universe expansion and in the centers of very massive stellar objects. 
Note that since we assume p > 0, then e > p and the equation of state for the 
second fluid ((5]) is below the stiffness limit. 

The gravitational equations for our two-fluid system 

R liV -\g» v R = k{Tl v + Tl v ) , (7) 

where k — 8irG/c is the Einstein constant, is equivalent to the Einstein equations 
in anti de Sitter space (with the negative cosmological constant —kp). Problems of 
finding exact wave solutions of Q within matter are in general very complicated, 
since one has to deal with interaction between the gravitational waves and the 
sound waves excited by th em. Even for the empty space-time there exists only few 
exact wave solutionjP 1 ^. As it was mentioned above, in this paper we consider 
the most studied wav e-like solution, which is based on the cylindrically symmetric 
Einstein-Rosen ansat^^^ 

ds 2 = e a( 7 -*) , dt 2 _ ^ _ r 2 e -W d6 2 _ e ^ dz 2 ^ (8) 

where the quantities 7 and ip are functions of the time t and the radial coordinate 
r only. 

As was shown irPH for any perfect fluid with a non-zero pressure, the energy- 
momentum tensor and its conservation can be written in terms of a scalar function. 
The non-rotating stiff fluid obeys particularly simple equations that arc equivalent 
to those of the massless scalar field a, 

' dvW^tf^dvo) = Oct = a - a" - -a' = , (9) 



s/-9 r 
where overdots mean time derivatives and primes stand for derivatives with respect 
of radial coordinate. 

Now we want to modify the representation used irPH for stiff matter to accom- 
modate it for our two fluid system. The simplest possibility is to shift e and p by p 
and describe the fluid characteristics by the scalar function a in the form: 

1 . ... iV'o 
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The right hand side of the Einstein equations ([7]) becomes equivalent to the energy- 
momentum tensor of scalar field, 

k (2u u u u - g U u)(p + p) = 2 d^odvO - g^ v d a ad a a . (11) 

So |(7j) may be written in the simple form: 

= 2d^ad v a . (12) 

For fluids that are not stiff the equations analogous to ^ and (fT2")) are more com- 
plicated and not amenable to the separation of variables employe d in th is paper. 
Non-zero components of the Ricci tensor for the ansatz ([8} ar jSSEQ] 

Rtt = $- i>" ~ i>'/r - 2^ 2 - 7 + 7" + i/r , 

R rr = -i> + ip" + ip'/r - 2i\)' 2 + 7 - 7" + i/r , 

R eg = -r 2 e- 2 ^ - 4>" - <///r) , (13) 

= e" 2 ^ - i>" - V//r) , 
i? rt = -2^/ + j/r . 
Then the Einstein equations in our case (fT!Zj) reduce to 

□a- = = , 

7 = 2r + , 

7' = r (^ 2 + V/ 2 + ct 2 + ff' 2 ) , (14) 

e -p = ?5 + p=I (a 2 - a' 2 ) e- 2 ^) . 

We see that both the fields, the fluid factor a and the metric component ip, 
behave like massless, minimally coupled, source-free scalar fields. Finding a and ip 
from the first equations of the system (fT4l) one can find expression of the second 
unknown metric parameter 7 in the formal 

7 = J dr r(& 2 + a' 2 + ij; 2 + 1P' 2 ) + / dt 2r{aa' + ipip') . (15) 



Some well studied special cases of the system ([14]) and (|T5|) arc: 
a). If we assume that the vacuum energy is equal to zero, p = 0, and also choose 
the solutions of (fl"4"|) in the form: 

a = at , ip = , (16) 

where a is a constant, we get singula r-fre e, static, cylindrically symmetric solution 
of Einstein's equations for stiff matteJ22l, 



ds 2 = e a2r2 (dt 2 - dr 2 ) - r 2 d9 2 - dz 
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b) . If again we suppose that p — and take 

(7 = 0, (18) 

the Einstein equations ([7]) reduce to another known syste m of e quations describing 
cylindrical gravitational waves in empty space (e 

= p= oj29130| 

$ - ip" - V = , 
r 

7 = 2ripip' , (19) 
7' = r (V 2 + V^ 2 ) • 

Now let us consider the case that we are looking for. Suppose the core of some 
supermassive object obeys the equation of state ((5]). The scalar modes of the stiff 
matter fluctuations a are indistinguishable from the excitations of metric ip. In this 
resonant case we write, 

V>(*,r) =a(t,r) , (20) 

and the main system (|14j) again reduces^ to the equations for cylindrical gravita- 
tional waves ()19|) . but waves are now inside our two-component medium (e,f>, p ^ 0). 
So we have cylindrical gravitational waves inside of almost stiff matter with the en- 
ergy density 

e = p + 2p = p + - vb a )e- 2 ^-^ . (21) 

Now we approaching the main task of our paper, finding stationary coupled 
matter-gravitational wave solutions of the system (fT4| , 



V>(r, t) = f(r) sm(wt) , a(r, t) = f(r) cos(ut) , (22) 
where, according to the first line of p^|) . f(r) satisfies the equation 

/" + -/' + w 2 / = . (23) 
r 

For (f22|) the second equation of (fl4|) gives the trivial relation: 

7 = 0. (24) 

Thus 7 = 7(r) is a function of the radial coordinate r only. Substituting ([22l into 
the third equation of (fT4|) we get: 



7' = r(c 2 / 2 + /' 2 ) . (25) 
Putting now (|22|) into the fourth equation of ([M]) one finds: 

e = p + 2p = p + \ [lo 2 sin 2 (c^t)/ 2 - cos 2 (u,t)f 2 ] e -^h(r)-^t)] t (26) 

where ip(r, t) and 7(7-) are defined via /(r) according to the relations (f2"2| and ([^5|l . 
respectively. 
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Note that the simplest solution of (f2"3")l is a zero order Bessel function of the first 
kind, 

f(r) = J (usr) , (27) 

which has the following asymptotics (see e.gPS): 

? ? 

US T 

r^O /r«l-^; 



r — > oo fir) w W cos fwr — — ) . (28) 

V 7rwr V 4/ 

Thus at the infinity [r — > oo), 

7 ^wr, e^p + e^, (29) 

and the energy e decays exponentially to the 'vacuum energy' value p. So (|22l) is 
a particula r case (z- independent solution) of Bessel waves extensively studied in 
recent year d 20 ! 21 ^. 

The coherent process described by the solution (12"2"1) m imics ver y much the be- 
haviour of electromagnetic waves in two-level medidjJ 2 ^^'^^'^^ , where the elec- 
tromagnetic waves are periodically absorbed and radiated by the two-level atoms, 
and thus the electromagnetic waves are trapped all the time inside the medium. 
Similarly, in our case gravitational waves are trapped by two-fluid system trans- 
ferring periodically the energy to the matter waves and then this energy restores 
again. 

According to the requirement that the energy in ([26| should be positive, it 
follows the restriction on the oscillation frequency 

us < n ~ ^/kp . (30) 

This itself means that the radial size of the object that traps nonlinear gravita- 
tional wave, as it follows from the same relation (|26p . should be larger than the 
characteristic length 

£1 \fkp~ ^ 
Supposing that kp is the ordinary cosmological constant we find unacceptably 
large the value of the p arameter ([3Tj) (of order of the horizon 10 28 cm), as astronom- 
ical observations implySSElEHl t ne present value of cosmological constant cannot 
exceed 10 -56 cm -2 . However, we introduce the fluid with the properties of the cosmo- 
logical constant as one of the ingredients of our two component system only for the 
simplicity of considerations. Vacuum energy is just needed as reference to keep posi- 
tive the energy of the system during coherent oscillations of gravitational and matter 
waves. Also problems of the definition of gravitational energy (in general) and equa- 
tions of states for very dense obj ects are w ell known. The controversial point about 
the concept of energy is the origin of a long-standing discussion 

on the energy and momentum carried by a gravitational wave. Perhaps one can find 
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other two (or more) component gravitating systems that will provide trapping of 
nonlinear gravitational waves of different kind by the mechanism considered in this 
article. 

Concluding we have investigated the coupled stationary state of a two-fluid 
gravitating system (stiff matter plus vacuum energy) with trapped Einstein-Rosen 
type nonlinear gravitational waves. In this medium the gravitational wave amplitude 
varies harmonically in time transferring the energy coherently to the stiff matter 
wave, and then the process goes to the backward direction, i.e. matter waves transfer 
the energy to the gravitational waves. We have defined the limits for the frequency 
of this energy transfer oscillations. 
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